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ABSTRACT 

We analyze the infrared (IR) spectral energy distributions (SEDs) for 10;um < A KSt < lOOyum for ~ 600 
galaxies at z ~ 0.7 in the extended Chandra Deep Field South by stacking their Spitzer 24, 70 and 160 /mi 
images. We place interesting constraints on the average IR SED shape in two bins: the brightest 25% of 
z ~ 0.7 galaxies detected at 24 pm, and the remaining 75% of individually-detected galaxies. Galaxies without 
individual detections at 24 pm were not well-detected at 70 pm and 160 pm even through stacking. We find that 
the average IR SEDs of z ~ 0.7 star-forming galaxies fall within the diversity of z ~ templates. While dust 
obscuration Llr/Luv seems to be only a function of star formation rate (SFR; ~ Lir + LuvX n °t of redshift, 
the dust temperature of star-forming galaxies (with SFR ~ lOM yr _I ) at a given IR luminosity was lower at 
z ~ 0.7 than today. We suggest an interpretation of this phenomenology in terms of dust geometry: intensely 
star-forming galaxies at z ~ are typically interacting, and host dense centrally-concentrated bursts of star 
formation and warm dust temperatures. At z ~ 0.7, the bulk of intensely star-forming galaxies are relatively 
undisturbed spirals and irregulars, and we postulate that they have large amounts of widespread lower-density 
star formation, yielding lower dust temperatures for a given IR luminosity. We recommend what IR SEDs are 
most suitable for modeling intermediate redshift galaxies with different SFRs. 
Subject headings: galaxies: evolution — galaxies: starburst — infrared: galaxies 



1. INTRODUCTION 

Dusty, intensely star-forming galaxies (SFR> lOM yr _1 ) 
are the dom inant contributors to the z > 0.5 cosmic SFR 
densit y (e.g.JFlores et al.lll999fc lElbaz et alJl2002t iPozzi et all 
I2004t iLe Floc'h et all 120051) . The thermal dust emission 
from these galaxies accounts for much of the cosmic IR 
background, which contains half of the radiation energy of 
the extragalactic background light ([Hauser & Dwek 1200 It : 
Lagac he et al.l l2005: Dole et al. 2006j). Therefore understand- 
ing the IR SEDs of the dusty star-forming galaxies is essential 
for mapping the evolution of star formation as a function of 
cosmic time, and is a key observational ingredient of our un- 
derstanding of galaxy evolution. 

Detailed studies of the IR SEDs over the full range of 
3-100 0yum have only been carried out for nearby galaxies 
(e.g.. lDaleet"al1l2005l) . finding t hat the IR SED shape is cor- 
related with dust temperature ((So ifer & Neugebauer Il99ll: 
Dale & Heloul 120021: Ichapman et all 120031: lLagache et all 
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20031) . Interestingly, despite this diversity of SED shapes, 
local galaxies' mid-IR (10-30/vm) luminosities are tightly 
correlated with their total IR lum i nosities with a scatter 
of ~ 0.3 dex (ICharv & Elbazl l200lt iPaoovich & Belli 12001 



iTakeuchi et alj2005btlDale et al.ll2005h 7 

This local observation has often been used as a key as- 
sumption in studies exploiting deep mid-IR imaging from the 
Infrared Space Observa tory at 15^m and the Spitzer Space 
Telescope at 24 pm (e.g JFlores et alJl999ll2004tlZheng eTaT] 
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lHammer et alJl2005t iBell et al.ll200H iMelbourne et all 
20051) . Such studies have found significantly enhanced star 
formation at 0.5 < z < 1, compared to the present day. 
Yet, this conclusion rests critically on the extent to which 
mid-IR luminosities reflect the total IR luminosity. Unfor- 
tunately, only a small fraction of mid-IR detected sources 
can be individually detected in the far-IR bands (e.g., Spitzer 
70 & 160//m), owing to limited signal-to-noise ratio (S/N ) 
and source confusion (iDole et alJ l2004b: Fra yer et alJ l2006). 
Thus, testing of this k ey assumption remains limited to the 
brightest sources (e. g., ISaiina et al.ll 2006 t Borvs et al.ll2006l) 
or is indirect (e.g.. lAppleton et al] | 2004t lYan et all 2005; 
iPope et alJl200alMarcillac et al.ll2006h . 

The goal of this paper is to explore the average IR SEDs of 
a stellar mass-limited sample z ~ 0.7 galaxies in the extended 
Chandra Deep Field-South (E-CDFS). In previous works, we 
have shown that stacking noise-limited images of a set of 
galaxies allows one to securely detect the mean flux of the 
galaxy set substantially below the individual detection limit 
dZheng et al.l2006l : lDole et al.l20 06). Here, we stack at longer 
wavelengths, i.e., 70 pm and 160 pm, to empirically deter- 
mine the population-averaged IR SED. We combine these re- 
sults with morphologies, and average fluxes at shorter wave- 
lengths, to place constraints on the dust extinction and SEDs 
of these galaxies. In §2 we describe the multi-wavelength data 
used to construct galaxy SEDs and the samples of z ~ 0.7 
galaxies. §3 presents our stacking methods. In particular 
we test the results of stacking noise and confusion limited 
70 and 160/vm images. In §4 we present the properties of 
the average SEDs. Discussion and conclusion are given in 
§5. Throughout this paper, we assume = 0.3, £2a = 0.7 and 
Hq =70kms _1 Mpc~' for a A-CDM cosmology. All magni- 
tudes are given in the Vega system except where otherwise 
specified. 

2. THE DATA AND SAMPLES 
2.1. The data 
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We use Spitzer 24, 70 and 160 pm data to study the thermal 
dust emission of z ~ 0.7 galaxies. In addition, we include 
deep ultra violet data from th e Galaxy Evolution Explorer 
(GALEX: iMartin et al]|2005al) . optical data from the Classi- 
fying Objects by Medium-Band Observations (COMBO-17; 
IWolf etal] 120031) survey, and four b and (3.6, 4.5, 5.8 and 
8.0 pm) Infrared Array Camera (IRAC; lFazio et alj2004 data 
to construct the stellar SED of a galaxy. These data cover 
wavelength range from 0.15 /mi to 160 pm in the observed 
frame, equal to the rest-frame range ~0.09 to ~ 100 pm for 
z = 0.7. 

The COMBO-17 survey has imaged the 30^5 x 30' E-CDFS 
in five broad (U, B, V, R and /) and 12 medium optical bands, 
providing high-quality astrometry (uncertainties ~ 071) based 
on a very deep /?-band image (26 mag at the 5<x li mit), photo- 
metri c redshifts (6z/(l + z ) ~ 0.02 at m R < 23; IWolf et alj 
120041) and stellar masses dBorch et all 120061) for ~ 11,000 
galaxies with ;«r < 24. We use the photometric redshift and 
stellar mass catalogs to select galaxy samples. 

GALEX ultraviolet observations provided deep far- 
ultraviolet (FUV; 1350-1750A) and near-ultraviolet (NUV; 
1750-2800A) images centered on the E-CDFS. The FUV and 
NUV images have a field of view of one square degree, a typi- 
cal point-spread function (PSF) of full width at half-maximum 
(FWHM) ~ 5", a resolution of 175 pixel -1 and a depth of 
3.63 pJy at the 5<x l evel. The data reducti on and source detec- 
tion is described in Morrissev et al.l d2005l) . 

The deep IRAC 3.6, 4.5, 5.8 and 8.0pm imaging data and 
MIPS 24, 70 and 160 pm imaging data were obtained as 
part o f the first run of MIPS GTO observations (Rieke et al. 
2004). A rectangular field of ~ 90' x 30' was observed in 
all bands (with small shifts between different bands). The 
effective exposure time is 500 s for the four IRAC band im- 
ages, 1378 spix -1 for MIPS 24 pm, 600 spix -1 for 70 pm and 
120spix -1 for 160 pm. IRAC 3.6 and 4.5 pm images have a 
PSF of FWHM ~ 178 and 5.8 and 8 .0 pm images have a PSF 
of FWHM ~ 270 dHuanget al.ll2004l) . The 24 pm image has a 
PSF of FWHM ^6". Sources ar e detected at 24 pm d own to 
83 p]y (at 80% completeness; see Papo vich et alj |2004 for de- 
tails of data reduction, source detection and photometry). The 
10pm image is characterized by a PSF of FWHM =* 18" and 
a resolution of 979 pixel -1 . The 160 pm image has a PSF of 
FWHM ^40" and a resolution of 16" pixel -1 . Sources with 
fluxes of fio >15 mJy can be individua lly resolved at 70 pm 
and of /i6o > 50mJy at 160 pm (see iDole etaill2004al for 
details). 

We take COMBO-17 astrometry as the reference coordinate 
and cross-correlate all other band catalogs with the COMBO- 
17 catalog. In each cross-correlation, we use bright stars and 
compact sources to estimate the systematic offsets and uncer- 
tainties between two coordinates. A position tolerance of 4cr 
uncertainty is adopted so that objects in the two catalogs hav- 
ing coincident coordinates within the tolerance (corrected for 
the systematic offsets) are identified as the same objects. The 
nearest COMBO-17 object is chosen if multiple ones exist 
within the tolerance. The adopted tolerances are 275 (FUV), 
370 (NUV), 170 (3.6pm), 172 (4.5pm), 175 (5.8 and 8.0pm) 
and 272 (24pm). For the 70 and 160pm catalogs, we firstly 
cross-correlate them with the 24 pm catalog as individually- 
detected 70 and 160 pm sources are bright at 24 pm; we then 
associate 70 and 160 pm objects with COMBO-17 objects us- 
ing the coordinates of the corresponding 24pm sources. The 
tolerance between 70/160 and 24pm is 5"/16" respectively. 
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Fig. 1. — The rest-frame color U - V versus stellar mass for a sample of 
1092 galaxies of known (GEMS) morphology in the redshift slice z = 0.7 ± 
0.05. Open symbols show the objects detected individually at 24 pm (/24 > 
83/jJy) and skeletal symbols represent individually-undetected ones. Most 
red (U - V >0.7) elliptical/lenticular galaxies and blue (U - V <0.7) small 
(M, <10M o ) spiral/irregular/peculiar galaxies are not individually detected 
at 24//m. The number of objects in each category is given in the brackets. 
The dashed line shows the mass cut M„ 
subsequent analysis. 
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2.2. The samples 

We combine spec troscopic redshifts fr om the VLT VI- 
MOS D eep Survey dLe Fevre et al. ll2005l) and the GOODS 
survey dVanzellaetal.ll2005L 12006) with photometr ic red- 
shifts from the COMBO-17 survey dWolf et alj 120041) to se- 
lect sample galaxies. The E-CDFS has HST ima ging from 
the G alaxy Evolution from Morphology and SEDs dRix et all 
2004, GEMS;) Survey. Visually classified morphologies are 
available for 1458 g alaxies with »ir < 24 in a thin redshift 
slice z = 0.7 + 0.05 dBell et alj|2005h . Of the 1458 galaxies, 
1 1 14 galaxies have all observations, including GALEX, IRAC 
and MIPS observations. To avoid contamination from active 
galactic nuclei (AGNs), we remove 22 X-ray detected s ources 
in the Chandra 250 ks observation (Leh mer et al.l l2005). leav- 
ing a sample of 1092 galaxies. Spectroscopic redshifts are 
available for 64% of the sample galaxies. The contribution 
from the X-ray-undetected AGNs to the total 24 luminos- 
ity of z < 1 galaxies is s uggested to be <~ 10% (Zheng et al. 
120061; [Brand et al.l2006l) . This will not have significant effects 
on our results. 

The sample demographics are shown in FigureQ] The sam- 
ple is limited by /?-band apparent magnitude (;«r < 24), cor- 
responding to approximately the rest-frame B-band at z ~ 0.7. 
Accordingly, the completeness of the sample, in terms of stel- 
lar mass, is a strong function of color: the mass limit for red 
(old or dusty) galaxies is M, ~ 10 10 M whereas blue galax- 
ies can be included down to almost M, ~ 10 9 M Q . We choose 
to impose a stellar mass cut of M* > 10 10 M G in what fol- 
lows; not only are almost all 24 pm emitters above this mass 
cut, but also the completeness of this sample is not a strong 
function of color (i.e., age or dust obscuration). We also se- 
lect another sample to explore the relationship between the 
IR SED shape and 24 pm luminosity; we extend the redshift 
slice to 0.6< z <0.8 to increase the number of sample galax- 



Infrared SEDs of z ~0J Star-Forming Galaxies 



3 





70 |j.m 


160 jam 

' I * 






4* 




» 

simulated 70 [J,m 

• * 

t 

" * 


Simulated 160 |im 

f 3 




• — - * • 

• 






* 






Fig. 2. — Observed MIPS 24, 70 and 160 /im image sections (15' X 15') in the extended CDFS and simulated 70 and 160 /jm images. The simulated 70 and 
1 60 jum images are created by degrading the 24 /im image to the 70 and 1 60 /.im resolution. No additional noise is added to the simulated images. 



ies. The final sample then comprises some 579 galaxies with 
AT. > 1O 1O M in the redshift range 0.6< z <0.8. Of them, 
218 are individually-detected at 24 /im with fluxes in excess 
of > 83/zJy; none is individually detected at 70 or 160 /im. 
X-ray detected sources have been excluded. 

3. IMAGE STACKING 

As outlined earlier, current missions are unable to yield 
individual detections for the vast majority of intermediate- 
redshift objects at far-IR wavelengths, owing both to contribu- 
tions from instrumental noise and confusion noise. In order to 
place constraints on the shape of the IR SEDs of 'typical' star- 
forming intermediate redshift galaxies, stacking on the posi- 
tions of known star-form ing galaxies can lower the effective 
noise (Zheng et al. 2006), allowing detection of the 'average' 
galaxy. In Zheng et al. (2006) we presented a description of 
stacking of 24 /im data (in that case to resolve the 24 /zm lu- 
minosity of dwarf galaxies); we briefly summarize the most 
important aspects of 24 /im stacking in 33.11 The focus of 
this paper is stacking at longer wavelengths, at 70 /im (PSF 
FWHM 18") and 160 /im (FWHM 40"), discussed in $H\ 

3.1. FUV, NUV and MIPS 24 /urn image stacking 

While much of this paper describes stacking results for sub- 
samples that are individually detected at 24 /im, some sub- 
samples are not individually detected at 24 //m. Further- 
more, the 24 /zm, FUV and NUV images share many of the 
same characteristics: the PSFs have similar FWHM, and at 
each wavelength, > 1/2 of the extragalactic background at 
that wavelength is resolved by these im ages (i.e., the images 
are only mildly confusion-limite d; e.g., Papovich et alj|2004l ; 
IXu et alj|2005t iDole et alj|2006l) . Accordingly, stacking of 



FU V, NUV and 24 /zm images is carried out in the same way 
(see lZheng et~ai1l2006l for more details). 

Three basic steps are adopted to derive the mean fluxes for 
galaxy subsets. First, we subtract all individually-detected 
sources from the ima ges. This is done u sing the software 
tool STARFINDER 7 jDiolaiti et all l2000h with an empiri- 
cal PSF constructed from 18/42/56 bright point sources at 
24 /zm/FU V/NUV respectively. Then we perform mean stack 
of the residual image postage stamps centered on the opti- 
cal coordinates of the objects that are individually undetected 
in the subset of interest. An aperture of 5" is used to in- 
tegrate the central flux of the mean-stacked image and esti- 
mate the background from the outer regions. Aperture correc- 
tions of a factor of 1.88/1.19/1.14, derived from the empiri- 
cal PSF at 24 /zm/FU V/NUV respectively, are adopted to cal- 
ibrate the stacked fluxes to the total fluxes. Last we sum the 
fluxes of individually-detected sources and the stack flux of 
individually-undetected sources in each galaxy subset, giving 
the mean 24 /zm/FU V/NUV fluxes. Uncertainties are derived 
from bootstrapping. 



3.2. MIPS 70 and 160 /zm image stacking 

The 70 /zm and 160/zm PSFs are considerably larger than 
those at shorter wavelengths, yielding confused images, re- 
solving only the brightest, relatively nearby sources — some 
<30% of the extr agalactic background at this wavelength 
dDole et al.ll2004al) . Galaxies at 0.6 < z < 0.8 are heav- 
ily confused in all but the brightest cases dDole et alj|2004at 
lLagache et aO2003ll2004l) . requiring stacking to gain insight 

7 STARFINDER gives identical results within the errors to the tool ALL- 
STAR in IRAF package. 
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Fig. 3. — Recovery of the mean fluxes of randomly distributed 24 /jm sources through stacking after convolution to the 70/im (left) and 160 /jm (right) resolution. 
The fluxes were taken from the 8255 individually-resolved sources in the ~ 1?5 X 0?5 24/im image and were split into 18 flux bins in the range from 83 /Jy to 
lOmJy. The number of objects in each bin decreases with increasing flux. Three measures of the stack flux are presented for each bin: the integrated central flux 
of the mean stack using the mean background for sky subtraction; the integrated central flux of the mean stack using the median background; and the integrated 
central flux of the median stack using the median of the background (see text for details). Errorbars are derived from bootstrapping. While the flux can be 
recovered with little bias at 70/;m resolution, flux recovery at 160/im leads to a systematic overestimate. 



TABLE 1 

The numbers of sources and average fluxes for 18 flux bins composed of 8255 individually-resolved 24/jm sources. 



No. 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 



</ 24 >(/iJy) 87 99 115 132 153 176 204 235 271 311 359 416 480 559 637 747 863 1946 
Number 1070 1259 1067 894 803 683 508 418 328 267 222 174 134 98 72 46 39 169 



into t heir long-wavelength IR SEDs (see, e.g., iDole et all 
12003) . 

The large angular extent of the long- wavelength PSFs poses 
a significant challenge for those wishing to estimate their av- 
erage properties. At z = 0.7, the PSF size of the 70 /im image 
(FWHM =i 18") corresponds to a physical scale of ~130kpc. 
For the 160 /mi image (PSF FWHM =; 40"), the correspond- 
ing scale is ~290kpc. Thus, the stacking results are a re- 
flection of the IR-luminosity weighted two-point correlation 
function on > lOOkpc scales, overestimating the true average 
fluxes of the galaxies of interest. 

In order to understand this source of systematic error in 
better detail, we carried out some simulations where syn- 
thetic 24 /mi data (using the observed positions and fluxes 
of individually-detected 24 /mi sources) are degraded in res- 
olution to the resolution of the 70 /im and 160 /mi data. For 
the purposes of this test, we assume a constant ratio between 
24 /mi flux and the wavelength of interest. We explore two 
cases. First, the positions are randomly scrambled (i.e., the 
relative brightnesses of galaxies are preserved but their posi- 
tions are random). This test gives an indication of the sys- 
tematic effects of stacking randomly-positioned sources. The 
second case is when both the fluxes and positions of sources 
are preserved. This second case is our 'best' estimate of the 
likely systematic uncertainties of stacking in a realistically- 
clustered case. Figure |2] shows the 24 /mi images degraded 
to the 70 and 160 /mi resolution respectively, compared with 
the observed MIPS 24, 70 and 160 //m images. Because of 
confusion, only a handful of bright 24 /im sources can be in- 
dividually resolved in the degraded images. Through stacking 
the mean flux can be estimated for 24 //m source subsets; the 



accuracy of the recovery of the mean flux shows how well 
stacking works at that corresponding image resolution. 

3.2.1. Stacking randomly distributed sources 

The whole 24 /mi mosaic image of the E-CDFS covers 
a rectangular sky area of ~ 1°5 x 0°5 and contains 8255 
individually-detected sources at 5cr detection limit (83 /iJy). 
The vast majority of the 8255 sources (fu > 83/iJy) are point 
sources. Replacing all sources with the 24 /mi PSF (empiri- 
cally constructed from bright stars), we generated an artificial 
24 fim image having the 8255 sources randomly distributed 
into a ~ 1°5 x Q°.5 blank field. The artificial 24/im image 
was then degraded to the 70 and 160//m image resolution 
respectively. The degraded images are then stacked at the 
24 /mi source positions, giving indications as to the expected 
quality of stacking results at 70/im and 160//m. We note 
that the replacement of the background blank field with the 
PSF-subtracted 24 /im image (having some contribution from 
individually-undetected sources) does not modify our results 
significantly in the mean, adding only some modest additional 
scatter. 

The 24/im sources were sorted int o 18 bins in flux , rang- 
ing from 83/iJy to lOmJy (see also IDole et alj|2006l) . with 
~ 1000 objects per bin at faint flux levels, and < 100 per bin 
at brighter fluxes. The numbers of sources and average 24 /im 
fluxes per bin are listed in Table Q] In real stacks, we choose 
to PSF-subtract out the individually-detected sources, stack- 
ing the remaining image (in order to reduce bias from bright 
sources in background estimates). In the real 70 and 160 /im 
images, we detect ~ 130 sources in each image. Thus, we 
must subtract approximately this many bright sources from 
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Fig. 4. — The recovery of the mean properties of 24/mi-detected sources at their observed positions. Left: Stacking results at 70/im resolution. Right. Stacking 
results at 160^m resolution. The errorbars reflect the variance from bootstrapping. 



the simulated images before stacking to avoid unrealistically 
biasing the results. Ideally, one would subtract all detected 
sources in the simulated 70 and 160 pm images; unfortu- 
nately, as the simulated and real 70 and 160 pm images have 
(unavoidably) somewhat different noise and brightness distri- 
bution properties, the number of detected sources is different, 
and substantially higher in the case of (the rather deeper) sim- 
ulated 70 pm data. Thus, we choose to PSF-subtract the 169 
brightest sources (the brightest of the 18 original bins) before 
stacking for the remaining 17 bins. 

For each flux bin, postage stamps were cut from the simu- 
lated image centered on the positions of the sources and were 
stacked. The size of postage stamps is 2'5x2'5 for the 70 pm 
image and 4' x4' for the 160/im image, allowing for a proper 
background estimate. Two stacked images were created for 
each bin through averaging or medianing. For the stacked 
70 pm postage stamp images, an aperture of radius 0'49 is 
used to integrate the central stack flux and an annulus with 
inner radius 0'82 and outer radius l'.23 to estimate the back- 
ground. For the stacked 160 pm postage stamp images, the 
corresponding aperture is of radius 1 '07 and the correspond- 
ing annulus is of inner radius If 07 and outer radius 1 '87. 
The mean and median of pixels in the annulus region were 
taken as the background of the mean-stacked image for sky 
subtraction. Only the median (nearly identical to the mean) 
was adopted as the background of the median-stacked image. 
Therefore, we derive three measures of the stack flux for each 
subsample: the integrated central flux of the mean stack using 
the mean background for sky subtraction; the integrated cen- 
tral flux of the mean stack using the median background; and 
the integrated central flux of the median stack using the me- 
dian of the background. Finally, aperture corrections derived 
from model PSFs (http://ssc.spitzer.caltech.edu/mips/psffits/) 
were applied to correct the estimates of stack flux to the total. 
The aperture correction is a factor of 1 .30 for 70 pm, 1 .25 for 
160 pm using the median background, and 1.34 for 160 pm 
using the mean background. 

Figure[3]shows the results of stacking randomly distributed 
sources at 70 and 160 pm resolution. Uncertainties are de- 
rived from bootstrapping. At 70 pm resolution, the mean flux 
of 24 pm sources of comparable flux can be properly recov- 
ered over a flux range of one order of magnitude. At 160pm 



resolution, recovered fluxes are of reasonable quality at bright 
limits, and become progressively more biased (by ~ 20-50%) 
towards fainter limits. After some investigation, it became 
clear that the the stack recovery at 160 pm is correlated with 
the number of objects in the stack bin. The increase of objects 
in number leads to an increase of the overlap between these 
objects within a fixed area. This increasing overlap leads to 
an increase of the stack flux relative to the input flux. By 
stacking a number of identical sources without background 
and foreground sources, we obtained similar results as shown 
in Figure [3] Indeed a point source contains 80% of its flux in 
an area of 3.6 arcmin 2 in the 160pm imaging and an area of 
0.9 arcmin 2 in the 70 pm imaging. To fill in a 1°5 x0°5 field, it 
requires around 750 160 pm sources or 3000 70 pm sources, 
compared to ~50 - 1260 objects in the stack bins. There- 
fore the stack results at 70 pm are little affected by the overlap 
between the stack objects but the stack results at 160pm are 
significantly influenced for stack bins of ~ 1000 objects. 

The three measures of the stack flux are nearly identical 
within the errorbars (~ 10%) for simulated 70 pm stacking and 
slightly different for simulated 160 pm stacking, in particular 
for the low flux bins, although the scatter is significant (~ 10 
- 25%). Note that we divided sources into stack bins in terms 
of their fluxes. For the stack bins having sources whose fluxes 
span a wide range, the median stack substantially underesti- 
mates the mean flux. 

3.2.2. Stacking sources at their observed positions 

In order to build a more realistic picture of the expected un- 
certainties from stacking at 70 pm and 160 pm, we repeat the 
last analysis except we keep both the fluxes and positions of 
the sources fixed to those observed at 24 pm. This gives in- 
sight into the influence of clustering on the results. In this test, 
the observed 24 pm map was degraded to the 70 and 160 pm 
resolution respectively, as shown in Figure [2] As previously, 
sources were split into 18 bins in 24 pm flux, and the stacking 
results for each flux bin were calculated. Figure |4] shows the 
results. Comparing the results with Figure [3] one can clearly 
see that the mean and median recovered fluxes are signifi- 
cantly affected by the clustering of sources on the sky. In 
general, our test suggests that stacking subsets of 70 pm and 
160 pm sources may overestimate the mean flux somewhat; 



6 



Zheng et al. 



10" - 



10 ,u ^ 



10 a - 



■ Irr/Pec (45) 
E/S0(17) 
Spiral (80) 



0.1 



1.0 10.0 
Rest Wavelength (urn) 



100.0 



Fig. 5. — The average SED, spanning a factor of 1000 in wavelength, as 
a function of morphological type for 24 /im detected (fi4 > 83/iJy) galaxies 
with M, > 10 10 M a in a thin redshift slice 0.65< z <0.75. Vertical errorbars 
are the lcr uncertainties derived from bootstrapping and horizonal errorbars 
show the band widths. The number of objects in each morphology bin is 
given in the brackets. 

typical biases are ~ 20%, with uncertainties of order ~ 20%. 
Yet, the magnitude and direction of these effects on different 
subsamples are difficult to estimate a priori. Figure[4]demon- 
strates that from bin to bin, the source clustering may lead to 
either an increase or a decrease of the mean flux of the source 
subsample, depending on the actual distribution and density 
of sources of different brightnesses on the sky. 

3.2.3. Stacking z ~0.7 galaxies 

The above tests demonstrate that it will be difficult to esti- 
mate the mean flux at 70 pm and 160 pm to better than ~ 20%, 
and that the extent of the bias towards higher or lower flux 
is difficult to predict with accuracy. In this paper, we adopt 
a pragmatic approach: for each subsample that is stacked at 
10 pm and 160 pun, we attempt to derive an individual 'bias 
correction' based on the 24 /mi image. Sources bright at 
24 pm are generally bright at 70 and 160/im, although the 
24 pm to 70 or 160/im flux rati o varies from obje ct to object 
at a scatter at the 0.5 dex level dDale et al.ll2.Q05b . Thus, we 
adopt the simplistic assumption that the 24 pm image can be 
taken as a good proxy for the 70 or 160 /mi image at ^ 6" res- 
olution and high S/N; 70 /mi and 160//m resolution images 
derived under this assumption (Fig. [2]) appear a reasonable 
description of the real 70 /mi and 160 /mi images. Thus, for 
a given subset of z ~ 0.7 galaxies, we estimated precisely 
the mean 24 /mi flux from their 24 //m images. Then we de- 
graded the 24 /mi images to 70 and 160 /mi image resolution. 
Two stack fluxes were estimated by stacking each of the two 
sets of degraded images. By comparing the stack fluxes to the 
actual mean 24 /mi flux of the galaxy subset, we obtained em- 
pirical corrections, which were applied to the corresponding 
70 and 160//m stack fluxes of the subset of z ~ 0.7 galax- 
ies, respectively. Uncertainties in this correction are applied 
also, in quadrature, to the derived stacking results at 70/mi 
and 160 /mi. 

4. RESULTS 

To explore the average IR SEDs of z ~ 0.7 star-forming 
galaxies, we first look into the dependence of the IR SED 
shape on galaxy morphology. Then, we investigate the rela- 
tionship between the IR SED shape and the 24 pm luminosity 
and the z ~ 0.7 IR SEDs to those of present-day star-forming 
galaxies. Finally, based on the IR SEDs, the extrapolation 
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Fig. 6. — Average SED as a function the 24 fim luminosity for galaxies with 
M* > 10 10 M G in the redshift range 0.6< z <0.8. The vertical errorbars show 
the lcr uncertainties and horizonal errorbars show the band widths. The 579 
galaxies are split into three bins in terms of their 24//m luminosities: L24- 
high, L24-medium and L24-I0W. The numbers in brackets denote the number 
of objects included in each bin. 

of the total IR luminosity from the 24 pm luminosity is dis- 
cussed. 

4.1. 77ie relationship between IR SED and morphology 

We used a sample of 1092 galaxies of known morphology 
to investigate the dependence of the IR SED shape on mor- 
phology. The 1092 galaxies are classified into three mor- 
phological types: elliptical/lenticular (E/S0), spiral, and ir- 
regular/peculiar. Figure Q] shows the relationship between the 
rest-frame color U — V and stellar mass. Objects detected 
at 24pm (J24 > 83/Jy) are marked with open symbols and 
undetected ones with skeletal symbols. The detection limit 
of 83 pJy corresponds to the monochromatic observed-frame 
24pm luminosity vL v> 24 = 6 x 10 10 L o at z = 0.7. The 24 pm- 
undetected galaxies are intrinsically faint in the IR bands (see 
34.21 >. and we excluded them in constructing these compos- 
ite SEDs. To avoid selection bias in color (see §12.21 i. we 
also excluded about one third of the 24 /mi detected galaxies 
with M* < 1O 1() M , leaving 152 galaxies in the final sample. 
We calculated the average luminosities in 14 bands for each 
subsample of 24 /mi-detected galaxies: FUV & NUV from 
GALEX, U, B, V, R and / from the COMBO- 17 survey, IRAC 
3.6, 4.5, 5.8 and 8.0/mi, MIPS 24, 70 and 160 pm bands. 
The average luminosity of each band includes contributions 
from both individually-detected and individually-undetected 
sources (at wavelengths other than 24pm). Errors in the av- 
erage luminosities were derived from bootstrapping, includ- 
ing contributions to the uncertainty from measurement errors 
in both the individually-detected fluxes and the individually- 
undetected fluxes. The average SED spans a range in the rest- 
frame of 0.1 to 100//m. 

Figure shows the average SEDs of the three subsam- 
ples. Each of the three average SEDs is dominated by a dust- 
extincted stellar spectrum at A < 5 pm and emission by dust 
at A > 5 pm. Irregular/peculiar galaxies show a higher ratio 
of dust to stellar emission compared to the spirals and E/S0 
galaxies. The slope of the IR SED (rest-frame 10-100/mi) 
is steeper for irregular/peculiar galaxies, somewhat interme- 
diate for E/S0 galaxies, and lowest for spirals although the 
uncertainties of 70 and 160/im luminosities are large. The ir- 
regular/peculiar galaxies usually form stars in relatively con- 
centrated regions, leading to a high star formation intensity 
(i.e. SFR per unit area). In contrast, the spirals are often char- 
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TABLE 2 

Average luminosities of massive-galaxy subsets in the redshift slice 0.6 < z < 0.8. 



Nobj 


log < M* > 
(Me) 


log vL v _24 
(Le) 


lo g yL v70 
(L e ) 


corr.(70) 


lo g vL v.l60 
(io) 


corr.(160) 


log < Ln— ioo > b 
(As) 


log < Lir > c 

(Le) 


58 

160 

361 


10.7 
10.5 
10.7 


10 52 +0 - 04 

10 09+ 001 
iu.u? 001 

Q T 1+0.03 

" -0.03 


10.79+°!° 

10 51 +016 
io.ji _ Q24 

o 90+0 28 
'■'"-0.99 


1.51 ±0.24 
1.03*0.32 
1.19±0.96 


II 17+O.I6 

III '-0.26 

10.98™-| 
< 10.27 


1.04±0.39 
0.92±0.61 


11 30 +01 ° 
-6.13 

11.04! | 


11 42+ 007 
-0.08 

11 12 +005 
10.28+» 2 | 



a The 70 fim and 160 //m luminosities estimated from stacking have included the corrections derived from degraded 24 fim images. The correc- 
tions are listed in the tabled L\2 ioo is the luminosity between 12 to 100 fim, calculated by linearly interpolating observed 24, 70 and 160 fjm 

luminosities (in logarithm space). Lir refers to the total IR luminosity between 8 to 1000/im. Local 1R SED templates from La gache et all 
(2004) Chary & Elbaz (2001) and Dale & Helou (2002) are used to fit the observed data points. The total IR luminosity is derived from the dust 
temperature-match templates. See text for details. 



acterized by a relatively low star formation intensity as the 
star-forming regions are widely distributed over disks. The 
star formation density for the E/S0 galaxies is somewhat be- 
tween those of the irregular/peculiar galaxies and spirals. The 
shapes of the IR SEDs are thus a function of star formation 
intensity in the sense that the dust temperature is primarily 
colder in systems of relatively lower star formation intensity 
(we will return this topic in ^5). 

4.2. The relationship between IR SED and 24 pm luminosity 

We used a mass-limited sample to study the relationship be- 
tween IR SED shape and 24 pm luminosity. This sample con- 
sists of 579 galaxies with M* > 10 10 M G and 0.6 < z <0.8. 
The sample galaxies were divided into three 24 pun lumi- 
nosity bins: Z.24-high, L24-medium and L24-I0W. The first 
two bins contain 218 individually-detected 24 pm sources 
(.724 >83yuJy) and all individually-undetected 24 pm sources 
(361 of the 579) are in the third bin. The Z.24-high bin is 
chosen to contain 58 brightest 24 fim sources so that its total 
24 fim luminosity equals that of the L24-medium bin. Con- 
sequently the stacked 70 and 160/im fluxes are expected to 
have comparable signal-to-noise ratios for the two bins. Aver- 
age luminosities in all 14 bands were calculated for the three 
subsets of galaxies. The average luminosities in the 24, 70 
and 160 fim bands are listed in Table [2] along with the num- 
ber of objects and mean stellar mass for each of the three sub- 
sets. The empirical corrections adopted for the 70 and 160 fim 
stack fluxes (see 33.2.31 1 are also presented. Errors include the 
uncertainties in measurements and bootstrapping errors. Fig- 
ure|6]shows the average SED from the rest-frame wavelengths 
0.1 to 100 /mi as a a function of the 24 fim luminosity. It is 
clear that the observed 24 fim (rest-frame 14 fim) luminosity 
is correlated with the 70 and 160//m IR luminosities for mas- 
sive galaxies at z ~ 0.7 in the sense that the 70 and 160 fim IR 
luminosities increase as the 24 fim luminosity increases; i.e., 
24 fim luminosity typically reflects high IR luminosity, rather 
than an enhanced rest- frame mid-IR excess (see also Bavouzet 
et al. 2007, in preperation). The Z.24-high bin has an average 
galaxy stellar mass 0.2 dex larger than the L24-medium bin; 
this can also be seen from the redder optical colors and higher 
rest-frame ~ 3 fim luminosity of the Z.24-high bin. The L24- 
low bin contains 361 massive galaxies that are individually 
undetected at 24 fim, including early-type galaxies with lit- 
tle star formation and late-type galaxies in the quiescent star 
formation phase (see Figure [TJ. The short-wavelength part 
of their average SED is dominated by a relatively old stellar 
population. The large errorbars of the 70 and 160/zm lumi- 
nosities compared to the small errorbar of the 24 fim luminos- 
ity are partially due to the intrinsic scatter among the sample 
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Fig. 7. — Relative dust temperature estimate from the IR flux ratio 
/60//100 versus f is/fm- Asterisks represent local star-forming galaxies from 
IDale etalJ <200Q) with IRAS observations at the 60 and 100pm and ISO 
observation at 15/im. The ratio fsal fioo is an indicator of the dust temper- 
ature, which tends to be correlated with IR luminosity in the local Universe. 
The dotted lines show typical values of feo/Zio o corresponding t he given 
log(Li R /L ) derived from IRAS data I Soifer & Neugebauer 1991). These 
lines are arbitrarily positioned on the Y-axis. The open circle shows the well- 
known local ultraluminous IR galaxy Arp 220. The squares show the three 
subsets of massive (M* > 10 10 M Q ) galaxies in redshift slice 0.6< z <0.8, 
labeled with the corresponding IR luminosities log(LiR/L ). Note that error- 
bars of the two flux ratios are not independent. The arrow indicates the effect 
of an increase in faa by a factor of 1.5. 

galaxies in this bin. 

4.3. Comparison with local SEDs 

A key goal of this paper is to compare the observed av- 
erage IR SED shapes at z ~ 0.7 to local 'template' SEDs. 
We adopted a s ample of local star-forming galaxies from 
iDale et al.1 d2000l) with IRAS and ISO 15 pm observations and 
total IR luminosity (8 - lOOO/vm) spanning from ~ 10 9 to 
10 12 Lq. Figure [7] shows these nearby galaxies in the IR flux 
ratio /6o //ioo versus f\sl fa) plot. The local galaxies distribute 
along a sequence with considerable scatter. The sequence 
is correlated with both dust temperature and IR luminosity. 
The dust temperat ure and IR luminosity incre ase for increas- 
ing /go/ /ioo (e.g.. ISoifer & Neuge bauerj[l991l) . The average 
IR SEDs for z ~ 0.7 galaxies are determined by the three 
MIPS bands, corresponding to the rest-frame ~14, ~41 and 
~94pm bands. The 24pm and 160/im luminosities can be 



8 



Zheng et al. 



1.5 



* 1.0 



0.5 



0.0 1 



10.0 E 





L 24 low (361) 



Lagache et al. (2004) 

Chary & Elbaz (2001) 

Dale & Helou (2002) 




-0.4 -0.3 -0.2 -0.1 0.0 0.1 -0.4 -0.3 



L 24 high (58) 



-0.2 -0.1 

log (f 6 (Aoo) 



> 



1.0: 



0.1 




0.0 0.1 -0.4 -0.3 -0.2 -0.1 0.0 0.1 



L 24 low (361) 



L04, log(L IR /y=11.4 
CE01, log^/L^1 1 .4 
DH02, log(L IR /L s )=11.3 



L04, log(L IR /L a )=11.1 
CE01, log(L IH /L s )=1 1 .0 
DH02, logfLm/L^II.I 



mo- 



1.0 r 



L04, log(L IR /y=10.2 

CE01, log(L IR /L5)=10.3 

DH02, log(L IR /y=10.3 




> 



0.1 




L04, log(f 60 /f 100 )= -0.31 

CE01, log(f 60 /f 100 )= -0.16 

DH02, log(f 60 /f 100 )= -0.28 



L04, log(f 50 /f 100 )= -0.36 

CE01,log(f 60 /f„, )=-0.38 

DH02, log(f 60 /f 100 )= -0.38 



L04, log(f 50 /f 100 )= -0.33 

- - CE01,log(f 60 /f 100 )=-0.41 
DH02, log(f 60 /f 100 )= -0.38 



10 



100 



10 100 
Rest Wavelength (urn) 




Fig. 8.— Fitting three sets of local templates from lLagache etalj 120041) . ICharv & Elbazl <2001l) and lDale & Helou 1 200^) to the average IR SEDs of three 
subsets of massive (M* > 10 10 M Q ) galaxies in the redshift slice 0.6< z <0.8. Only the three MIPS data points (24, 70 and 160/im) of each observed IR SED 
are used in the fitting. Top: The panels show^- 2 as a function of the template's characteristic log(/6o//ioo)- For each observed IR SED two best-fit templates are 
chosen: one is matched in luminosity (i.e., with the normalization constant equal to unity) and the other is dust temperature match (or SED shape match; i.e., 
with the minimum^ 2 ). Bottom: The upper panels shows the best-fit luminosity match templates and the lower panels show the best-fit dust temperature match 
templates, compared to the observed SEDs from the rest-frame 2 to 100 pm delineated by IRAC and MIPS data points for the L24-high bin (left panels), the 
Z.24 -medium bin (middle panels) and the L24-I0W bin (right panels). The SEDs are normalized to unity at 2.1 pm. Note that the stellar components of the locoal 
templates are somewhat arbitrarily set. The disagreements between the averaged SEDs and the best-fit templates at A rcst < 5pm should be ignored. 



taken as rest-frame 15 fim and 100//m luminosities, for which 
'^-corrections' are negligible. We estimated the rest-frame 
60 fim luminosity by linear interpolation between the MIPS 
measurements in log-log space (as shown in Figure |6]). We 
used the local sample to test the linear interpolation. First, we 
derive the 41 fim fluxes by linear interpolation between IRAS 
25 and 60 /mi measurements in log-log space. Second, we de- 
rive 60 fim fluxes by the same method between 41 and 100 fim 
for the 59 local galaxies adopted. The estimated 60 fim fluxes 
are 25+1 1% lower than the observed IRAS 60 jum fluxes. This 
hints that our rest-frame 60 fim luminosities of the z ~ 0.7 
galaxies might be underestimated. We compared the z ~ 0.7 
galaxies to the local galaxies. As shown in Figure UJ the 
two populations are roughly located in the same region of the 
/60//100 versus /15//60 plane. Specifically, the z ~ 0.7 star- 
forming galaxies of IR luminosity 11 < log(LiR/L ) < 11.4 



(those in the L24-high bin and L24-medium bin; the estimates 
of the total IR luminosities will be discussed later) populate 
the relatively low-temperature end of the template sequence. 
In the local Universe, these low-temperature galaxies tend to 
be of relatively low luminosity log(LiR/L Q ) < 10.5. This, with 
significant uncertainties, suggests that the typical dust temper- 
ature of z ~ 0.7 luminous IR galaxies (LIRGs; i.e., galaxies 
with log(LiR/L ) > 1 1) is lower than that of local galaxies of 
comparable IR luminosity. The IR flux ratios are poorly de- 
termined for the L24-I0W bin because of the large uncertainties 
in the average 70 and 160//m IR luminosities, which are par- 
tially due to the intrinsic scatter among the sample galaxies in 
the L24-I0W bin (including early-type galaxies and late-type 
galaxies in the quiescent star formation phase). 

We compare d our average S EDs with the IR SED mode l 
templates from lLagache et al.1 d2004l) . ICharv & Elbazl d2001l) 
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Fig. 9. — The ratio between the rest-frame 15 fim to the 12-100 /im luminosity (left plot) and to the total IR luminosity (8-1000 /im; right plot) as a function of 
the IR flux ratio /so//loo- Asterisks are local galaxies collected from literature with IRAS, ISO 15 fim observations and at least one observation at wavelengths 
longer than 100 /an. Squares show our results for galaxies with M* > 10 10 M Q and 0.6< z <0.8. The relations derived from SED templates are presented for 
comparison. One can see that the extrapolation of the rest-frame 15 pm luminosity based on the three sets of templates give comparable estimates of total IR 
luminosity within a considerable scatter. Our data points distribute within the scatter of the local star-forming galaxies but towards the colder templates. Roughly 
speaking, a factor of 10 is a quite good average bolometric correction to the estimate of the total IR luminosity from the rest-frame 15 fim luminosity. Generally 
speaking, the use of local templates of Ljr < 10 10 - 5 L Q for extrapolation of total IR luminosity from intermediate redshift 24 fim luminosities gives a more 
accurate result than the use of Lir ~ 10" L Q templates. 



and iDale & Heloul (120021) . The three sets of templates were 
empirically calibrated to represent local star-forming galaxies 
spanning a wide range in the IR flux ratio feo/fioo- Dale & 
Helou's IR SED templates are charact erized by the flux ra- 
tio fml fioa- We use the equation from Chap man et al.1 (120031) 
to parameterize the IR luminosity (8 - lOOO/um) as a func- 
tion of /6o//ioo- 8 The tem plates from Lag ache et al.l (120041) 
and ICharv & Elbazl (200 1 ) were characterized by the IR lu- 
minosity. We extended each set of SED templates by linear 
interpolation in logarithmic space to a grid of SEDs with the 
characteristic IR luminosity ranging from 10 9 to 10 13 L o with 
a resolution of 0. 1 dex. The observed IR SED of each subset 
was compared to the three sets of templates: 



X 



z 



Lobsj f scale X Lt&mp,i\Z — 0.7) 



(i) 



where L b s .;, L temp ,; and cr, are the observed and template lumi- 
nosities and their uncertainty in filter i, respectively, and / sca ] e 
is a normalization constant. Here only 24, 70 and 160/im 
bands were used to fit templates, i.e., Affiiters = 3. L t emp is cal- 
culated by convolving the redshifted SED template to z = 0.7 
with the 24, 70 or 160 /im filter transmission function. / sca ] e is 
chosen to minimize the^- 2 for each template. The top plot in 
Figure [8] shows the fitting results. The templates with a nor- 
malization constant / sca i e of unity are chosen as the best-fit lu- 
minosity match templates and those with a minimum^- 2 as the 
best-fit dust temperature match (or SED shape match) tem- 
plates. The best-fit templates are compared to the observed 
SEDs in the bottom plot of Figure [8] 
As shown in Figure [8] the local SED templates track the 



8 The total IR luminosity defined in Chapman et al. 12003) is between 3 
and 1 100 , in good ag reement with the adopted one IDale et al.l200lt sec 
also|Takeuchi et al 20053). 



observed SEDs reasonably well, in particular for the L24- 
medium bin, which contains the majority of intense star- 
forming galaxies. As the total IR luminosity is dominated 
by emission in the rest-frame wavelength range between 10 
and 100//m, the estimates of the total IR luminosity with tem- 
plates from different models give similar results. In contrast, 
X 2 is a measure of the shape agreement between a template 
and an observed SED. It is worth noting that the estimate of 
total IR luminosity based on the three measurements at 24, 
70 and 160/im is not sensitive to the shapes of the IR SED 
templates. Both the best-fit luminosity match and the best-fit 
dust temperature match the SED templates and suggest a total 
IR (8 - 1000 fim) luminosity of log L 1R /L Q -11.4, 11.1, 10.3 
for the three average SEDs, respectively. Generally speak- 
ing, the shape (or dust temperature) of the average IR SEDs 
of star-forming galaxies at z ~ 0.7 (i.e., the L24-high bin and 
L24-medium bin) are better fitted by the local SED templates 
of characteristic IR luminosity Lir < 10" L than those of 
Lir > 10 n L . This holds for all three sets of templates. It 
confirms that the typical star-forming galaxies at z ~ 0.7 are 
likely to have relatively colder dust emission than local galax- 
ies with comparable IR luminosity. 

4.4. The extrapolation from the rest-frame 15 fim to the total 
IR luminosity 

Local IR SED templates are often used to estimate the to- 
tal IR luminosities from single mid-IR band luminosities for 
distant star-forming galaxies. With the averaged IR SEDs de- 
termined in the three MIPS bands, we are able to better con- 
strain the estimates of the IR luminosities. We compared the 
estimates with those transformed from single 15 fim luminosi- 
ties using local SED templates. 

By linearly integrating between the average 24, 70 and 
160/im luminosities (in logarithm space) for the star-forming 
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Fig. 10. — The ratio IR/UV flux (i.e., the fractional dust obscuration) as a 
function of SFR. Solid circles show the three subsets of galaxies with M. > 
10 10 M & and 0.6 < z < 0.8, sorted by the 24 /im luminosity. Other symbols 
show three morphology subsets of star-forming (f'24 > 83,ujy) galaxies in 
redshift slice 0.65 < z < 0.75. The dashed line shows the local relation 
I Mart in et alJ2005bl) . 

galaxies at z ~0.7, we estimated the rest-frame 12-100 /zm 
luminosities. Uncertainties are calculated from the combi- 
nation of the uncertainties in the three bands. The total IR 
(the rest-frame 8-1000 /mi) luminosity is derived from the 
three dust temperature-match templates to each observed IR 
SED (the 24, 70 and 160 /mi). We combine the scatter be- 
tween the templates and the uncertainty in determining the 
dust temperature-match templates, in which the uncertainties 
of the 24, 70 and 160 /mi luminosities are counted, as the un- 
certainty for the total IR luminosity. The results are listed 
in Table [2] For comparison, we use the NAS A/IPAC Extra- 
galactic database (NED) to collect a sample of 29 local star- 
forming galaxies (z < 0.1) with observations at 12, 25, 60 and 
100/zm by IRAS, at 15 /mi by ISO and at least in one band 
longer than 100 /mi. The 12-100 /mi luminosity and the total 
IR luminosity are calculated as above. When observations do 
not reach 1000 /mi, a modest (< 20% of the IR luminosity) 
extrapolation is employed using the local SED templates. 

Figure [9] shows the extrapolations (bolometric corrections) 
from the rest-frame 15 /mi to the 12-100 /mi luminosity (left 
plot) and to the total IR luminosity (right plot) as functions 
of the IR flux ratio /eo/Zioo- The local star-f orming galax- 
ies exhibit a sign ificant scatter dCharv & Elbaz 2001; see also 
iDale et al] [2005): such scatter is typically adopted as the sys- 
tematic uncertainty of template-based estimates of total IR 
luminosity from rest-frame mid-IR luminosity. The z ~ 0.7 
star-forming galaxies distribute within the scatter of the lo- 
cal star-forming galaxies. Again, there is a tendency for the 
z ~ 0.7 galaxies to cluster towards the colder templates, sug- 
gesting that use of spiral galaxy templates for extrapolation 
of total IR luminosity from intermediate redshift 24 /zm fluxes 
gives a more accurate result than the use of starburst/LIRG 
templates. 

4.5. Dust extinction 

A universal relation between SFR a nd dust extinction is 
suggested by several studies at z < 1 ( iHopkins et alj I20011 : 



l Adelbereer & Steidelll2lMlBel]||2003k l Takeuchi et alJ[2 005a: 
Zheng et alJl2006l:lBuat et al.ll2006t although see lReddy et al] 
2006 for a dissenting view at z ~ 2). With a measured esti- 
mate of the total IR luminosity from the 24, 70 and 160 /zm lu- 
minosities, we explore the relationship between SFR and dust 
extinction at z ~ 0.7, comparing it to that at z ~ 0. The dust 
extinction is described by the IR to UV ratio Lir/Luv- The 
UV luminosity, i.e., the integrated luminosity between rest- 
frame 1500-2800 A, is estimated from linear interpolation of 
the FUV, NUV, U, and B band luminosities fo r z ~ 0.7 
galaxy subsets shown in Figures[5]and[6] Following Be ll et al] 
(2005), we derived the SFR from the IR and UV luminosities 
with the formula 

SFR/(M G yr 1 ) = 9.8 x 10" n (Li R + 2.2L UV ), (2) 
assuming a stellar population with a constant SFR for 
100 Myr and a Rroupa initial mass function. Figure [10] 
shows the relationship between dust obscuration and SFR for 
z ~ 0.7 galaxies, compared to the local relation derived from 
IRAS and GALEX data with a sca t ter at ~0.5-ldex level 
(iMartin et al] 1200511 IXueFaT] l2006t iBuat et alj l2006h . As 
shown in Figure \W[ the z ~ 0.7 galaxies distribute perfectly 
along the local relation over one order of magnitude. 

5. DISCUSSION AND CONCLUSION 

We selected a sample of 152 galaxies in the redshift slice 

0. 65 < z < 0.75 of known morphology from HST imag- 
ing and another sample of 579 mass- limited (M, > 10 10 M G ) 
galaxies in the redshift slice 0.6 < z < 0.8 to study the IR 
SEDs at that redshift. We divided our sample galaxies into dif- 
ferent mass-limited morphology and 24 //m luminosity bins. 
For each bin, we determined the average luminosities in 14 
bands from the FUV to the FIR by summing the individual 
detections and adding in the stacked flux from non-detections. 
Careful efforts were taken in stacking the noise and confusion 
limited 70 and 160 /zm images. Empirical corrections, deter- 
mined from the 24 /zm image, were introduced to account for 
the clustering effects on the stack results. 

The average luminosities in three MIPS bands determine 
the IR SED from the rest-frame 10 to 100/zm. Our principal 
result is that the average IR SED shape of z ~ 0.7 intensely 
star-forming galaxies (with IR luminosities ~ 10 n L G ) is simi- 
lar to reasonably 'cool' local templates (i.e., templates of 'nor- 
mal' spiral galaxies). The dust SED seems to depend on mor- 
phological type for star-forming galaxies at z ~ 0.7. This has 
the immediate and important implication that the use of local 
templates to extrapolate total IR luminosity from observed- 
frame 24 /zm data is a well-posed problem, at least, on aver- 
age. Interestingly, galaxies with 'cool' dust temperatures in 
the local Universe all tend to have IR luminosities < 1O 1O 5 L , 

1. e., distant intensely star-forming galaxies tend to be charac- 
terized by colder dust emission than their local counterparts 
of comparable IR luminosity. 

Previous studies have found evidence for a somewhat 
cooler dust SED at 0.2 < z < 2.5 than for local galaxies of 
a com parable luminosity (e.g.. iPope et al.ll2006t ISajina et al] 
2006) at 0.2 < z < 2.5. These studies were selected in rest- 
frame > 100 /zm emission, and the authors suspected that their 
overall tendency towards 'colder' IR SEDs was in part due 
to that long-wavelength selection. Our sample is selected on 
rest-frame ~ 15 /zm emission, i.e., by warm dust; yet, we find 
a 'colder' average SED at a given luminosity than is found 
locally. This tends to support the interpretation that the off- 
set which we and others have found towards colder tempera- 
tures at a given luminosity are at least in part a real difference. 
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Fig. 1 1 . — Same as Fig. [7] The filled circles show the two subsets of galaxies 
in the redshift slice 0.6 < z < 0.8 split by 24 luminosity surface density. 
The inner panel shows the average SEDs of the two subsets. The XL24-high 
bin and the El24-1ow bin data points are located separated along the dust- 
temperature sequence of local star-forming galaxies, indicating that the £l24- 
high bin is characterized by a hotter dust emission than the El24-1°w bin. 
This is clearly seen from the difference between their IR SEDs. 

We suggest that this tendency toward colder dust tempera- 
ture reflects a difference in dust and star formation geometry: 
whereas local LIRGs tend to be interacting systems with rela- 
tively compa ct very intense sta r formation and of comparable 
masses (e.g. JWang et al.ll2 006). 0.5 < z < 1 LIR Gs tend to be 
disk-d ominated, relativ ely undisturbed g alaxies (IZheng et al.1 
12004 iBell etaT]|2005t iMelbourne et al.ll2005l) . We suggest 
that these disk galaxies host widespread intense star forma- 
tion, like star formation in local spirals but scaled up, lead- 
ing to rela tively cold du s t tem peratures. This is consistent 
with what Cha nial et alj (120061) found — the scatter in the 
^ir _ Tdust relation for star-forming galaxies is largely in- 
duced by the size dispersion of the star-forming regions; and 
the IR luminosity surface density Sir - r^st relation is more 
fundamental than the the Lir - 7d ust relation. We examined 
the IR flux ratio /60//100 (indicator of dust temperature) as 
a function of the IR luminosity surface density for z ~ 0.7 
star-forming galaxies. By taking gal axy size (i.e. half-li ght 
radius) derived from HST imaging (Hau61ereLafl[2007) as 
a proxy of the size of star-forming regions in a galaxy and 
24 jjm luminosity as a proxy of total IR luminosity, we split 



a sample of 3 1 1 24 /mi-detected galaxies in the redshift slice 
0.6 < z < 0.8 into two bins in the 24-fim luminosity surface 
density Sl24- The average SEDs for the two subsets of galax- 
ies were constructed in the same way as described in 34. H and 
shown in the inner panel of Figure [TT] The total IR luminos- 
ity is estimated as log(LiR/L ) =11.3 and 11.0 for the £l24- 
high bin and the Xl24-1ow bin respectively. As shown in Fig- 
ure [TT] the two subpopulations are distributed along the dust- 
temperature sequence of local star-forming galaxies. Galax- 
ies in the EL24-high bin on average show hotter dust emission 
than those in the Xl24-1ow bin, suggesting that the IR luminos- 
ity surface density p lays an essential role in shaping IR SED 
(IChanial et al.ll2006h . 

Finally, with UV luminosities derived from GALEX data 
and IR luminosities derived from three band MIPS data, we 
determined a quite precise relationship between the dust ob- 
scuration and SFR for z ~0.7. Our results show an excel- 
lent agreement between the SFR-dust obscuration relation at 
z ~0.7 with that at the present day, indicating that no signif- 
icant evolution occurs since that redshift. Our measurements 
give the mean values of SFR and dust obscuration. The ac- 
tual scatter in dust obscuration can spread by 1 - 2 dex for in- 
dividual galaxies of given SFRdMartrn et al.il2005bt IXu et all 
I2006t iBuat et al.ll2006b . iReddv et al.l d2006l) claimed that the 
dust obscuration for star-forming galaxies is systematically 
smaller at z ~ 2 than the present day. Their sample is dom- 
inated by rest-frame UV-selected galaxies. This may lead to 
a poten tial selectio n bias against objects with high IR/UV ra- 
tio (see IBuat et a l. 2006 for the comparison between UV and 
FIR selected samples of local galaxies and a similar discus- 
sion). Therefore studies based on unbiased samples of high- 
redshift galaxies will help to answer the question whether 
the SFR-dust obscuration relation still holds at z > 1 . This 
will add important constraints to our understanding of galaxy 
evolution involvi ng star formation and metallicity enrichment 
(IZhenget al.1120061) . 
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